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ABSTRACT 

Transverse magnetohydrodynamic (MHD) waves have been shown to be nbiq- 
nitons in the solar atmosphere and can in principle carry snfficient energy to gen¬ 
erate and maintain the Snn’s million-degree onter atmosphere or corona. How¬ 
ever, direct evidence of the dissipation process of these waves and snbseqnent 
heating has not yet been directly observed. Here we report on high spatial, tem¬ 
poral, and spectral resolntion observations of a solar prominence that show a 
compelling signatnre of so-called resonant absorption, a long hypothesized mech¬ 
anism to efficiently convert and dissipate transverse wave energy into heat. Aside 
from coherence in the transverse direction, onr observations show telltale phase 
differences aronnd 180° between transverse motions in the plane-of-sky and line- 
of-sight velocities of the oscillating hne strnctnres or threads, and also snggest 
signihcant heating from chromospheric to higher temperatnres. Comparison with 
advanced nnmerical simnlations snpport a scenario in which transverse oscilla¬ 
tions trigger a Kelvin-Helmholtz instability (KHI) at the bonndaries of oscillating 
threads via resonant absorption. This instability leads to nnmerons thin cnrrent 
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sheets in which wave energy is dissipated and plasma is heated. Our results 
provide direct evidence for wave-related heating in action, one of the candidate 
coronal heating mechanisms. 

Subject headings: waves — Sun: chromosphere — Sun: transition region 


1. Introduction 


Since 15 years, transverse oscill ations have been d etected quite regularly in the solar 


corona. The TRACE observations ( 

Handv e 

al. 

19991) starting from 1998 showed plentv 

of flare-induced oscillations ( 

Nakariakov et al. 

1999 

: Schriiver et al. 

1999; 

Aschwanden et al. 


19991) . The y were successfully used for seismology of the solar coro na, by estimating the mag - 
netic held flNakariakov fc Ofmanll200ll). the loon ra dial structure flAschwanden et al.ll2003ll. 


the v ertical density scale height flAndries et al.ll2005f) and Alfven transit times flArregui et ah 
20071) . 


Since 2007, with the advent of the CoMP instrument and the Hinode satellite flKosugi et ah 


20071 ). we are now convinced that the e ntire solar atmosphere is hlled with low-arnplitude 


propa gating and standing Alfvenic waves flTomczvk et al.ll2007l:l0kamoto et ah 


2007; 


De Pontieu et a 


2007h . even in the absence of a hare excitation. Subsequent studies with CoMP flTomczvk fc Mclntoshl 
2 OO 9 I) . Hinode and the Solar Dynamics Observatory {SDO; Pesnell et al.ll2012) spacecrafts re - 
vealed th e wave properti e s and showed their presence in prominences fISchmieder et al.1120131). 


spicules flHe et a 


201ll:IWang et al 


12009 


2012 


Okamoto fc De Pontieu 2011) and corona 


Nisticb et al 


2013 


Anhnogentov et al.1120131). These oscillations are 


the manifestation of magnetohydrodynamic (MHD) waves fe.g.. lOliverl 120091: iMackav et al. 


loops flMcIntosh et al. 


2 OIOI: lArregui et ahl 120121: iMathioudakis et al.l l2013l. and references therein). Such waves 


may play an important role in the so-called coronal heati n g problem by supplying th e 


Sun’s hot outer atmosphere with energy (e.g., lAlfvenI 119471: iParnell fc De Moortell 120121) . 


corona and accelerate the solar wind ( 

Uchida & 

<!abi] 

rakil 1974: Antolin & Shibata 2010: 

van Ballegooiien et al. 2011: Matsumoto & Suzuki 

2014 

). However, although such waves are 

sometimes damped on short spatial and temporal scales in observations ( 

Aschwanden et al. 


19991: iPascoe et al.ll2010L 120111) . it is still unclear whether any signihcant dissipation (and 
thus heating) occurs in the corona, which has very low resistivity. 

The dissipat ion of transverse MHD waves in an inhomogeneous medium has been studied 
for decades (e.g., Ilonson Il978l: iHollweg et al.lll990l: ISakurai et al.lll99ll: iGoossens et al.lll992 


Van Doorsselaere et al.l 120041) . Theoretically it is expected that the transverse, coherent 
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motion is subject to d amping by resona nt absorption in an inhomogeneous prominence thread 

Resonant absorption converts the tra nsverse wave 


(Arreeui et ah 

2008; 

Soler et ah 

2012) 


into azimuthal motions with a hne spatial scale in the tube’s boundary layer fjVerth et ah 


F K 

2010l:lSoler et al.ll2010[ lArregui et alJl2011h . The resonant and amplihed azimuthal motions 
are prone to be unstable to the Kelvin-Helmholtz instability (KHI), because of the large 
shear motions. The KHI deforms the boundary layer and leads to enhanced dissipation of 


the wave energy into heat in thin, turbulent current s heets flKarpen et ahlll993l:l0fman et ah 


1994 iLapenta fc Knoll! 120031: iTerradas et ahN2008a|) . Recent simulations bv lAntolin et ah 


1 ) 2014 ') show that even small amplitude oscillations can lead to the instability, whose vortices 
combined with line-of-sight (LOS) effects result in strand-like structures within coronal loops, 
but at scales that are difficult to resolve with current coronal instruments. 

This study is composed of two parts. In t his hrst paper, unique coordinated obser¬ 


vations with the Sola r Optical Telescope (SOT; iTsuneta et al.l I 2 OO 8 I: IShimizu et al.l 12008 


Suematsu et al.ll200^ of Hinode and NAS A’s recently launched Interface Region Imaging 
Spectrograph (IRIS: \De Pontieu et al.ll2014j) satellite are used to provide evidence for a tell¬ 
tale sign of resonant absorption and associated heating i n an active reg i on sol ar prominence. 
We extend the three-dimensional MHD model used in lAntolin et al.l (120141 ) to a scenario 
matching our observations. Through radiative transfer, we show that the observed trans¬ 
verse dynamical coherence, the specihc phase relations between the azimuthal motions and 
displacement of the threads and the characteristic heating are a characteristic feature of 
resonant absorption combined with the KHI. The details of the numerical model and a more 
thorough analy sis of the obse r vation al signatures of transverse MHD waves in prominences 
are provided in lAntolin et al.l (120151) . hereafter Paper 2. 


Observation and Data Reduction 


IRIS observed an active region prominence on the southeast limb of the Sun from 09:00 
to 11:00 UT on 19 October 2013 (Figured]). Four-step sparse raster scans, which have 1" 
gap between neighboring slit positions, were performed together with the slit-jaw images of 
the Si IV and Mg II lines. The cadence was 10 s for the slit-jaw images and 20 s for spectra 
at each position with the spatial and spectral resolutions of 0.33"-0.40" (240-290 km on the 
Sun) and 3 km s“^. The exposure time was 4 s in both. Coordinated observations with the 
SOT were performed from 09:09 to 10:27 UT. The SOT obtained Ca II H-line hltergrams 
with a cadence of 8 s, exposures of 1.2 s, and the spatial r esolution of 0.2" (150 km). The 
Atmospheric Imaging Assembly (AIA; iLemen et al.l l2012l) aboard SDO provided context 
hltergrams of multiple wavelengths with a cadence of 12 s. 
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Here we used level-2 data of IRIS, which are calibrated by dark current subtraction, 
flat-held correction, and geometrical correction. Remaining ohsets by a few pixels between 
the slit-jaw images and spectral data were found, and then we shifted the images to align 
with hducial marks that appear on the slit-jaw images and the spectra. Next, we per¬ 
formed co-alignment of time series of the slit-jaw images. Large-amplitude orbital variation 
of IRIS pointings, whic h is called wobble, is typically removed at the moment of planning 
(iDe Pontieu et al.N2014i) . but small misalignment less than 1” still remained. With cross 
correlation of two neighboring-time images of Si IV, the ohset information was derived and 
applied for all IRIS slit-jaw and spectral data. Hence, one can hnd slow drifts of the slit 
locations in the time-series of images (see the online movie A). Although this impact was 
small enough in a timescale of focused oscillations (10-20 min), we considered it when we 
chose the nearest locations to investigate the LOS velocity of oscillating threads. Lastly, we 
performed co-alignments between the IRIS slit-jaw images (Si IV) and the SOT filtergrams 
(Ca II) by adjusting their plate scales, rotating, and shifting images. To emphasize the faint 
prominence structures, a radial-density filter was applied for the disk parts (also part of 
spicules) of all images. No gamma value is applied for these images. 

We mainly use Mg II k (2796 A, typical of plasma at 10,000 K) spectra of the prominence 
to determine their LOS velocity as a function of time. Even though the spectra show some 
central reversal, we perform c entroiding to determin e the LOS velocity which provides a 
good estimate of the velocity (ILeenaarts et aLN2013h . Single Gaussian functions are used 
to £t the Mg II spectra to determine the centroid LOS velocity. Comparison of these LOS 
velocity signals with oscillatory displacements of the same threads in the plane-of-sky (POS), 
obtained in the Ca II H (3968 A, 10,000 K) passband of the SOT, allows to reconstruct the 
three-dimensional flow generated by the waves. 


3. Methods and Results 

3.1. Height dependence of threads and signatures of heating 

The Ca II H-line movie shows that the prominence consists of numerous horizontal 
threads that move horizontally with various speeds (10-40 km s“^). Space-time plots tracking 
the flow reveal that many of these threads oscillate in the POS, with signatures of damping. 
These features can also be seen in the Mg II k and Si IV (1403 A, 80,000K) IRIS images, 
despite the coarser spatial resolution. In the space-time plot of LOS velocity derived from 
Mg II k spectra along one slit location, we can find oscillatory patterns at all heights (Figure 
121(b)). Moreover, the upper region shows higher velocity amplitudes than the lower region. 
Similarly, the line width of Mg II k is larger in the upper region than in the lower region 
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(Figure [2]^c)). From the histogram of LOS velocity (Figure [2](d)), the higher region (regions 
3 and 4, indicated by yellow and red lines, respectively) has more high-speed components 
with more than 15 km s“^, while velocity in the lower region is more conhned around 0 km 
s“^. This is consistent with the large-amplitude motions seen in the Ca II movie. Similarly, 
the histogram of line width (Figure [2](e)) also shows stronger broadening in the higher region 
than in the lower region. We note that the line width in Figure |2] includes the instrument line 
broadening (5.5 km s“^ in FWHM). The movie shows that the lower region has more LOS 
superposition with a multitude of threads overlapping. One would expect such superposition 
to lead to broader spectral line prohles. However, this is not found, suggesting that the 
difference in non-thermal line broadening between the low and high regions is instead caused 
by increased turbulence or unresolved wave motions, which can lead to heating. 

Further careful investigation reveals that the horizontal threads in Ca II H have smaller 
length and shorter lifetime at higher altitudes. As the moving threads in the Ca II line fade 
away, co-spatial threads appear in the hotter Si IV line with similar horizontal speeds. Space- 
time plots of the Ca II and Si IV images show this more clearly (Figure [3]). Even though 
some locations have coexistence of cooler and hotter materials, trajectories of several threads 
show the transition from lower to higher temperature. These fade out hnally. This indicates 
that the cool plasma is heated to coronal temperature through the mid-temperature range 
of the Si IV line. 

Two kinds of intensity variability can be observed in the online movie A. First is a 
variability produced by the dynamic nature of threads linked to flows and waves. This 
effect, as discussed in the next section, is mostly behind the fact that a hxed slit can only 
capture parts of the transverse oscillations in the PCS, and entails a timescale on the order of 
one period or less. The second kind of variability corresponds to a gradual intensity change 
from the chromospheric emission in Ca II H to emission in the transition region range of the 
Si IV line, and can be appreciated in Figure [21 This variability is therefore linked to heating 
and occurs on slightly longer timescales of 1 — 2 periods. Both effects, dynamics and heating, 
are intertwined and make the exact determination of lifetimes of threads difficult, even more 
from the strong line-of-sight projection effects (especially in Ca II H). In the online movie 
A we observe threads appearing and disappearing in the top part of the prominence on 
timescales of 10 — 15 min, i.e. a couple of periods of the oscillations. Such lifetimes appear 
to be longer in the lower part of the prominence, where projection effects are stronger but 
also where dynamics are reduced, as explained previously. 
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3.2. Transverse oscillations 


The Ca II movie clearly shows the ubiquity of transverse motions of threads. Oscilla¬ 
tions over several periods are not easy to spot in this kind of dynamic prominence. This 
i s largely due to the horizontal flow, which makes the oscillations appear non-stationary 
flOkamoto et ahl 1200711 . and move them quickly away from the hxed slit positions. Also, as 
shown by Figure [Hand the Ca II movie, vertical displacements of large horizontal sections are 
observed. Both effects complicate the visibility of mnltiple periods. To clearly indicate the 
presence of transverse oscillations and better visnalize the longer periodic oscillations of the 
threads, space-time plots at different horizontal locations following the flow in the threads 
are provided in the online movies C-F. One can clearly hnd 2- to 3-period oscillations of 
PI and P3 threads with damping of their amplitndes (marked by arrows in the movies). 
Althongh the abnndance of oscillations is also clear for P2 and P4, the large quantity of 
threads and their very dynamic natnre impedes a fully clean tracking of oscillations beyond 
2 periods. Signatnres of damping are detected, but also constant and even some cases of 
increasing amplitnde oscillations. 


3.3. Characteristic phase differences and dynamical coherence 


Next, we focus on four oscillating cool threads that can be seen clearly and individually 
with little superposition in the Ca II images (P1-P4 in Figure [T](a)) to investigate the phase 
relationship between the transverse displacement in the Ca II line and LOS velocity in the 
Mg II line (Figure H]). We show the combined velocity signal from the nearest slit location 
at each time. Despite the fact that the threads oscillate across the IRIS slit locations at all 
times, the LOS velocity signals detected around the threads are strongly coherent (Figure 
ED. Interestingly, the oscillating threads show phase differences aronnd 90° and 180° between 
the transverse displacement in the POS and the LOS velocity. Large phase differences above 
90°, in particular 180°, do not match theories of the classic kink mode with discontinuous 


boundary in the c 
waves or rotation 


ensity prohle. nor th ose invoking classic axisymmetric torsional Alfven 


Goossens et ahl 120141) . 


Figure Elshows the locations of the slits aronnd the threads for positions P1-P4. Besides 
the ont-of-phase relation between the POS motion of the thread and the LOS velocity, an¬ 
other striking featnre is the existence of coherence in the Doppler signal across a signihcant 
transverse distance on the order of a thread diameter or more, away (above and below) the 
Ca II H thread. Such coherence may come as a surprise if one assnmes that the external 
medinm dominates our signal - if that were the case one would expect a loss of coherence. 
Besides coherence in the LOS velocity, the existence of coherent transverse POS motion can 
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also be clearly seen in the online movies C-F. Snch coherent POS motion of t hread-like strnc- 


tnre seems to be a common characteristic of transv erse waves in prominences fjOkamoto et ah 


o 

o 

CM 

Lin et ah 

2009; 

Lin 

2011 ; 

Hillier et ah 

2013 

) and in loops, as t 

(Ofman & Want 

>; 2008fl or coronal rain ( 

Antolin & Verwichte 

2011 ) 


tions we show evidence for the hrst time that the coherent motion extends to the azimnthal 
flow generated by the transverse wave. 


It is important to notice that the Ca II H threads appear snrronnded by Mg II k emission, 
as evidenced by Figure [U This fact, together with the existence of dynamical coherence in 
the POS and along the LOS at signihcant distance across the threads suggests the existence 
of larger prominence flux tubes containing the threads. This scenario will be further explored 
in the next section. 


3.4. The transverse MHD wave model 


To understand the observational results and in particular the peculiar phase difference, 
we use three-dimensional MHD simulations of a transversely oscillating prominence flux tube, 
combined with radiative transfer modeling. The flux tube has initial internal to external 
electron number density (nj/ue) and temperature (Tj/Te) ratios of 10 and 1/100, respectively. 
The transition layer between the internal and external medium has a width l/R ^ 0.4, 
where R is the flux tube radius. For the forwar d modeling we give typical values found 


in pr ominences n* = 10^° cm ^ and T* = 10^ K flTandberg-Hanssenl I1995 k IVial fc Engvold 


20151 1. At time f = 0 the flux tube is subject to a transverse perturbation generating a 
fundamental standing transverse MHD wave with longitudinal wavenumber k ps 0.015/i?, 
amplitude Vq =8 km s“^ and phase speed of 776 km s“^. The external plasma-is 0.01 and 
the fully ionized external medium has a magnetic held value of 18.6 G. We assume a long 
(L = 200 Mm) and thin (R = 1 Mm) prominence hux tube, matching the observations. For 
further details please see Paper 2. 

As expected from theory, resonant absorption in the simulation sets in immediately and 
the energy from the purely transverse kink mode is rapidly transferred to the azimuthal 
waves in the resonant layer at the boundary of the tube (the online movie B). This induced 
azimuthal component of the velocity is, at hrst, in phase with the dipole-like azimuthal how 
outside the tube, thu s exhibiting a 90° phase diherence with the transverse displacement 
flGoossens et al.ll2014j ). Since the Alfven speed increases as the density decreases the further 
away we are from the tube’s axis, but especially across the boundary layer, the periods of 
oscillation of the azimuthal waves in the resonant layer will be smaller than that of the 
transverse displacement, leading to a drift in time between the transverse displacement and 









































LOS velocity (a process known as phase mixing). Over long time scales of several periods 
this process leads to generation of small-scales from which a broadening of the line width, a 
decrease of the Doppler signal and a loss of coherence is obtained across the resonant layer. 
In the short observable window of a few periods following an external perturbation (before 
resonant absorption dampens out the transverse POS motion) a favored 180° phase difference 
stands out. This phenomenon is more clearly seen with a standing wave, for which the phase 
difference cannot grow indehnitely due to the limited length of the flux tube (Figure [6l see 
also the schematic representation in Figure [71(a)). 


Apart from resonant absorption, another important m echanism in our simulations is the 
Kelv in-Helmholtz instability (KHI). As in the coronal case ( Terradas et ahlboOSa : Antolin et al. 


20141) . the KHI is rapidly triggered at the boundaries after one period of oscillation, leading 
to the generation of vortices and current sheets (see Figure 2 in Paper 2). Importantly, a 
complex interplay between resonant absorption and KHI sets in, in which the instability 
extracts energy from the resonant layer and imparts momentum on the generated eddies, 
which rapidly degenerate into turbulent-like flows, as can be seen in the online movie B. 


To analyze the overall effect of resonant absorption combined with KHI we performed 
detailed radiative transfer modeling, explained in the next section. 


3.5. Radiative transfer 


The chromospheric conditions of our flux tube demand a completely differe nt approach 
in ter ms of forward modeling with respect to the previously modeled coronal case (jAntolin et al 
20141 ). In order to properly compare with observations we translate the numerical results 
into observable quantities by means of appropriate radiative transfer modeling. For opti¬ 
cally thin transition re gion lines such as Si IV 1402.77 A w e calculate the synthetic emission 
using the FoMo co de (lA ntolin fc V an Do orsselaerel l2013l) . based on the CHIANTI atomic 
database version 7 f Pere et al.f]l997n2009 ). For optically thick lines such as Mg II h&k and 
Ca II Hfc K we use the RH c ode, as follows. We employed the t wo-dimensional version of th e 
RH code (jUitenbroekI l200l[ based on the method described in iRvbicki &: Hummer! (1199211 ) 
to calculate the emergent spectra of the Mg II h&k as well as the Ca II H&K resonance lines 
in non-LTE through a slice perpendicular to the a xis of the prominence. Th e RH code uses 
the efficient Accelerated Lambda Iteration scheme flRvbicki fc Hummerlll992[) for transitions 
with overlapping wavelengths. It uses a formal solut ion of the transfer equation based on th e 
short-characteristics method flKunasz fc Auerlll9881) with Bezier interpolation flAuerll2003[) . 


Since the MHD code does not provide explicit electron number densities (required to 
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calculate the collisional excitation, de-excitation and ionization-recombination rates), they 
were computed with the RH code in a two-step process. First, an estimate of electron 
densities was made under the assumption of LTE. Second, this hrst estimate was used to 
calculate a full two-dimensional non-LTE solution of a six-level plus continuum hydrogen 
atom. Finally, the ionization of hydrogen (i.e., the proton density) gives an estimate of the 
non-LTE ionization contribution of hydrogen, which was added to the LTE contribution of 
all other elements (including He). 


We used abundances from iGrevesse fc AndersI fjl991l l for our calculations, which cor¬ 
respond to photospheric v a.lues. For the at omic models we employed the hve-level plus 
continuum model for Ca II fjUitenbroeklll989f) . the ten-level plus continuum model for Mg II 
f Uitenbroeklll997l). an d a hve-level plus continuum model for Hydrogen with transition values 


fjVernazza et al.lll98ll) . treating all pertinent radiative bound-bound and bound-free transi¬ 
tions explicitly. Since the prominence cross section is optically thick only in the very cores of 
the lines, we used the approximation of complete frequency redistribution (CRD) in all line 
transitions, including the Lyman lines in hydrogen. This is justihed because wavelengths 
outside the Doppler core are optically thin and scattering in them does not signihcantly 
affect the line transfer in the line cores. 


3.6. RH results 

Through forward modeling of the numerical results we generate observable quantities for 
direct comparison with the IRIS and Hinode observations. We hrst calculate the synthetic 
emission in Mg II h&k and Ca II H&K with RH on a two-dimensional cross-section at the 
center of the tube at a given time during the numerical simulation (we pick time t = 1 ,162 
s, corresponding to roughly two full periods from the start of the simulation; the numerical 
model at this time is shown in Figure 2 in Paper 2). In Figure [8]we show the source function 
at the center wavelengths of the Mg II k and Ca II H lines, as well as the contours for 
diherent values of the optical thickness. As can be seen, the source function shows a ring 
shape with maximum values around the boundary layer of the tube, and negligible values at 
the center of the tube. This means that the emergent intensity in both of these lines comes 
from the boundary of such flux tubes. 
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4. Discussion 

4.1. A ring shape for the source function 

The obtained ring shape for the source function in Mg II h&k and Ca II H&K is mostly 
a temperature effect. Cool prominence plasma around the boundary layer is heated up by 
the KHI to transition region temperatures, thereby enlarging inwardly the transitional layer 
from the core to the surrounding corona. The higher temperature in the outer ring raises 
both the electron density and the electron temperature, both giving rise to more collisional 
excitations in the h & k lines (that are then followed by an increased rate of spontaneous 
emission). The radiative transfer results further suggest that cool chromospheric cores of 
such flux tubes do not contribute to intensity in Mg II h&k but only to opacity. The opacity 
in these lines depends on the ground level populations in Mg II and not on the radiation 
field. Also, the emerging intensity is optically thin and is therefore mostly dependent on 
temperature. Accordingly, rays going through the center of the tube are optically thick (in 
k more than in h, by a factor of 2), so these suffer from self absorption of the signal that 
comes from the back surface. Such rays have an intensity ratio of k and h lines tending to 
1.5. On the other hand, rays going through the surface of the tube only have ratios close 
to 2 because they are purely optically thin. Ratios between 1.5 and 2 are thus obtained for 
spatially broad rays (simulating the case of coarse instrument spatial resolution). 

An important result from our forward modeling with RH is that the emergent intensity 
in both the Mg II and Ca II H lines comes from the boundary of such flux tubes, where 
resonant absorption and the KHI vortices take place. Furthermore, the optical thickness 
values in Figure [8] show that the boundary of such flux tubes is optically thin in both lines, 
while the cool and dense core is significantly opaque only in Mg II h&k, an important 
difference with respect to Ca II H&K. This result implies that the former lines suffer less 
from line-of-sight projection effects than the latter, and are therefore more favorable for 
prominence analysis at the limb such as that performed here. 

This result allows us to perform an approximation to optically thick radiation with 
FoMo, as explained in Paper 2, which retains the key aspects of the radiative transfer cal¬ 
culations. This approximation allows us to forward model swiftly and thus study in detail 
the observational signatures of our numerical model. 
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4.2. Interpretation of the observational results: resonant absorption &; KHI 


While the phase drift produced by resonant absorption that results in the transverse 
MHD wave model (which would correspond to the initial part of phase mixing) could poten¬ 
tially explain the ont-of-phase difference found in onr observations, it also poses a problem 
of detectability, and especially it fails to explain the strong dynamical coherence explained 
in section 13.31 Indeed, resonant absorption is expected to be conhned to the fine resonant 
layer, which is completely invisible to the cnrrent instrumentation. Even in inhomogeneous 
flux tubes with irregular internal density profiles, while the system still evolves with a main 
oscillating freqnency and a large-scale resonance pattern sets in around the ensemble of 
inhom ogeneities, the complicated resonant layer s are still fonnd to be highly spatially con¬ 
fined f Pascoe et ah 2011 : Terradas et ah 2008bl) (this case is sketched in Figure [7] (b) and 
we refer to it as the mnltiple kink wave scenario). Snch systems could potentially explain 
the POS dynamical coherency (assnming that the inhomogeneities within flux tubes could 
correspond to the observed threads) but would fail at explaining the LOS velocity coherency 
and especially the out-of-phase behavior. 


This is where the second mechanism in onr simulations, the Kelvin-Helmholtz instability 
(KHI) becomes important. Indeed, as explained in section [3~4I the instability extracts energy 
from the resonant layer and imparts momentnm on the generated eddies, which rapidly 
degenerate into tnrbulent-like flows that are still embedded in the large-scale azimnthal 
flow. This results in a spatial broadening of the transition layer between the prominence 
core and the external environment, which mostly retains the characteristic phase drift of 
the resonance layer with respect to the transverse displacement of the prominence core. 
Furthermore, the layer affected by the KHI, which occupies most of the original flux tube, 
is subject to significant heating, leading to a fade out of Mg II k emission, and gradnal 
appearance of Si IV emission. This is shown in Figure |6l where the Doppler signal captnred 
by the slits (placed a distance of a diameter away from the flux tube, above and below) shows 
the characteristic ont-of-phase behavior with respect to the POS motion of the prominence 
core. As explained in Paper 2, this is not only the result of the KHI broadening, but also 
of the PSF width of IRIS (which has been taken into acconnt in onr forward modelling), 
which spreads the signal over a significant distance. Furthermore, due to the significantly 
large perpendicnlar distance over which the LOS velocity is retained and to the symmetric 
Doppler profile with respect to the flux tube axis, the dynamical coherence in the observations 
is recovered. 


As shown in Paper 2, the characteristic observational signatnres of the transverse MHD 
wave model can be seen for basically any LOS in a plane perpendicular to the tube axis, 
and for optically thin and moderately thick conditions (Paper 2). In the mnltiple kink wave 
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scenario previously discussed the same observational features would be expected since the 
ensemble of ‘threads’ (enhanced density regions within the large flux tube) evolve coherently, 
developing a large azimuthal flow around the system. The onset of the KHI around this flux 
tube would therefore also produce the same characteristic features as in the present case. 


Another feature predicted by our numerical model is the appearance of thread-like 
structure in both chromospheric and transition region lines, which result from the KHI 
vortices and line-of-sight projection effects. Such structure can b e seen in Figure E] a nd is 
similar to the strand-like structure obtained in the coronal case (lAntolin et al.N2014l l. In 
Paper 2 we discuss the possibility of the observed Ca II H threads being the KHI vortices. 


4.3. Other interpretations 

Our detailed modeling shows that the observations are fully compatible with a state- 
of-the-art numerical model of transverse MHD oscillations and thus our interpretation of 
resonant absorption. Furthermore, we can also prove that in the scenario of other possible 
processes, such as rotation and the torsional Alfven wave, the observed coherence and phase 
difference between the POS motion and LOS velocity cannot be reproduced. 

Let us hrst consider the case of a torsional Alfven wave. We can further divide this 
case into 2 different cases that would lead to the observed POS oscillation of threads. Case 

1 (sketched in Figure [71(c)) corresponds to a scenario in which the ensemble of threads is 
subject to an m = 0 torsional Alfven wave (as belonging to a larger flux tube). Case 2 
(sketched in Figure [71(d)) shows a single thread and corresponds to a scenario in which each 
thread is subject to a coupled torsional Alfven wave and kink wave, in which the azimuthal 
flow from the former dominates that of the latter one. In each sketch a generic location for 
the IRIS slits has been taken, in which the width of the slit is set according to the size of 
the thread, as is observed. Case I can be further divided into two cases, depending on the 
size of the flux tube with respect to the size of the region scanned by the rasters. If the 
slit positions cover different parts of the flux tube (edges and center, as pictured by case I 
within Figure [71(c)) the LOS velocity signals recovered by the green and red slit positions 
should be out of phase, disagreeing with the observed transverse coherence in the Doppler 
signals shown in Figure [5l (in each panel of this hgure the signal from a slit follows the same 
behavior regardless of whether it is placed above or below the thread). If the slit positions 
cover only one half of the flux tube the transverse coherence is recovered (as shown by case 

2 in panel (c)), as well as the 180 degree phase difference. However, in this case a gradual 
increase of amplitude should be observed between the slits the closer they are to the flux 
tube boundary, which is not observed in our Hinode-IRIS measurements, as shown in Figure 
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|5] in all panels. For the case of conpled torsional Alfven and kink waves the Doppler shift 
will drastically change sign depending on where the slit is across the thread (as in the hrst 
case of the torsional Alfven wave), which impedes the collective in-phase signals between 
the slits that we observe in Fign reElin all panel s (asym metric \-y profile, cf. section 4.3 in 
Paper 2). This is also shown bv iGoossens et ahl fl2014j) . in their Fignre 9. 


Let us now consider the case of rotation. Similarly, we can identify two cases leading 
to the observed POS oscillation. The first case (sketched in Figure [71(e)) corresponds to 
a scenario in which the ensemble of threads is rotating coherently. In the second case 
(sketched in Figure [71(f)) each thread is subject to rotation coupled to a kink wave, in which 
the azimuthal flow of the former dominates over that of the latter one. In both cases a slit 
placed always on one side of the thread should not change sign in time. Also, a slit that 
crosses the thread should exhibit a change in the sign of the Doppler signal. None of the 
observed cases in Figure 5 satisfies both of these scenarios, which invalidates rotation as well. 


We note that the previous predictions should not be affected by optical thickness. In¬ 
deed, for both cases of rotation and a torsional Alfven wave, the frontside and backside of 
the flux tube with respect to the observer exhibit a coherent motion. Therefore, regardless 
of how opaque the core of the flux tube is, the signal detected by IRIS will retain its Doppler 
shift value. The predictions should therefore hold, excluding both scenarios of rotation and 
torsional Alfven wave. 


Figure [91 summarizes the observed features together with the expected observational 
features of the 6 considered cases. The plus/minus sign denotes a change in the LOS velocity. 
The green and red curves denote the shape of the LOS velocity in time, as detected by the 
slits of same color in the sketches. 


5. Conclusions 


Observations o f coronal transverse osci llations were hrst (and are often) reported in the 
aftermath of hares (iNakariakov et al.[[l999h . Such reports have been the subject of several 
dozen follow-up theoretical and observational studies, which suggest r esonant absorption as 
the rnost plausible mechanism causin g the damping of such oscillations (jPe Moortel fc Nakariakov 
2012c [Nakariakov fc Verwichte[[2005h . However, to the date no solid consensus exists on the 
causes of this damping. Furthermore these large-scale oscillations occur only in the aftermath 
of an impulsive event such as a hare. Such events cannot be the main mechanism of heating 
in the quiescent corona as they are too infrequent (compared to the typical coronal cooling 
times). Heating from Alfvenic waves invoking resonant absorption in prominences or coronal 
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loops has also prev iously been reported flNakariakov et al.lll999l: lOfman fc Mouradianlll996 


Ofman et al.lll998l) . supported by numerous theoretical studies suggesting t hat this mech¬ 


anism is efficient and ubiquitous flArregui et al.ll2008l: iGoossens et al.lI 2 OIII: iTerradas et al. 


2008bl) . However, such observational reports have been largely based on indirect evidence 


due to incomplete measurements (lack of spectroscopic information allowing the detection of 
the three-dimensional flow and heating generated by the wave) or lack of resolution leading to 
a precarious causal relation between resonant absorption and the subsequent heating. In the 
present work we show for the hrst time direct evidence of resonant absorption as measured 
with a spectrograph and associated heating to transition region temperatures. This asso¬ 
ciation is performed with unprecedented detail and therefore narrows down to unparalleled 
levels the causal relation in evidence for wave heating. 


Our numerical modeling suggests that an important coupling between resonant absorp¬ 
tion and KHI takes place, in which KHI efficiently extracts energy from the resonant layer 
and converts it into heat. As the KHI sets in viscous and ohmic heating stron gly increase due 


to th e generation of vortices and current sheets around most of the flux tube (IGoossens et al. 


2 OO 2 I ). As explained previously, the heating occurs rapidly at hrst and proceeds to a roughly 


constant rate for most of the simulation. The boundary layer, characterized by transition 
region temperatures around 10^ K, expands rapidly in the timescale of one period, mostly 
inwards towards the center of the tube, which leads to fading with time of threads in Ga II 
and Mg II (both 10,000 K), and the subsequent appearance (in the same location) of threads 
at higher temperatures on a timescale of 10-15 min (e.g.. Si IV at 80,000 K) (Figure [H]). 
Such heating is also present in our observations on similar timescales: as the moving threads 
in the Ga II line fade away (the online movie A), co-spatial threads appear in the hotter 
Si IV line with similar horizontal speeds. Space-time plots of the Ga II and Si IV images 
show this clearly (Figure [3]). Even though some locations show coexistence of cooler and 
hotter materials, the trajectories of several threads show a transition from lower to higher 
temperatures with threads hnally fading out even from Si IV images. This scenario cannot 
be explained by a multi-thermal structure of the threads, and is strongly suggestive of cool 
plasma being heated to coronal temperatures through the mid-temperature range of the Si IV 
line. Evidence for coronal counterparts of the prominence comes from the online movie A: 
while clear hne structures cannot be identihed in the Fe IX (~1 MK) channel of FDO/AIA, 
horizontally-elongated bright features do occupy the same region as in the cooler passbands 
(Figure [T]). 


The observed heating is directly predicted by our transverse MHD wave model and thus 
provides supporting evidence for our interpretation, which, as mentioned, is focused on the 
spatial coherence and the phase relationships between LOS velocity and POS oscillations. 
The combination of observations and modeling thus provide compelling evidence for wave 
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dissipation through the combination of resonant absorption and KHI and resulting heating 
to at least transition region temperatures in a prominence. The observed oscillations have an 
energy flux of 100-1,000 W m“^, where the velocity amplitude is 20 km s“^ in the case that 

the plasma density is typically 10“^^ 


the v iew angle is 45 


flLabrosse et ah 


10 kg m ^ 

2010l ). and the Alfven speed is assumed to be 300-600 km s“^. Here we have considered the 


more general case of the multiple kink wave scenario (sketched in Figure [71(b)), for which 
the energy fl ux is within a factor o f 2 from that of the bulk Al fven wave multiplied by a 
hlling factor flGoossens et al.ll2013l: IVan Doorsselaere et ahll2014f) . We take a conservative 
Filing factor of 5% to avoid overe stimation. The energy fl ux is enough to locally heat the 
plasma to coronal temperatures fjWithbroe fc NovesHlOTTl) . Our results may also pr ovide 
insight into the puzzling rapid disappearance of prominences fjTandberg-Hanssenlll995l) and 
i n the heating mec hanism of hot plasma bubbles that have been reported in prominences 
flBerger et al.ll201ll) . More generally, our results provide a pathhnder to search for heating 
resulting from resonant absorption in other solar regions (e.g., quiet corona and coronal 
holes) that are permeated with transverse magnetic waves. 


Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
domestic partner and NASA and STFC (UK) as international partners. It is operated by 
these agencies in cooperation with ESA and NSC (Norway). IRIS is a NASA small explorer 
mission developed and operated by LMSAL with mission operations executed at NASA 
Ames Research center and major contributions to downlink communications funded by NSC 
through an ESA PRODEX contract. SDO is part of NASA’s Living With a Star Program. 
Numerical computations were carried out on Cray XC30 at the Center for Computational 
Astrophysics, NAOJ. T.J.O thanks Ippon-kakou-kai for their encouragement and was sup¬ 
ported by JSPS KAKENHI Crant Number 25800120 (PI; T.J.O.). T.J.O and P.A. were 
supported by JSPS KAKENHI Grant Number 25220703 (PI: S. Tsuneta). B.D.P. was sup¬ 
ported by NASA under contract NNG09FA40C (IRIS), NNX11AN98G and NNM12AB40P. 
T.V.D. was supported by FWO Vlaanderen’s Odysseus programme, GOA-2015-014 (KU 
Leuven) and the lAP P7/08 CHARM (Belspo). 
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(a) CallH (c) Mg II k 



Fig. 1.— Prominence observed on the southeast solar limb on 19 October 2013. (a) Hin- 
ode/SOT image of the observed prominence in the Ca II H line (~10,000 K). The horizontal 
elongated structures are prominence threads. The vertical ones are spicules, which are not 
studied here. The four white lines are the positions of the IRIS slit. The green box indi¬ 
cates the region used for tracking the horizontally-moving threads for Figure [HI The vertical 
red bars are the selected locations for the space-time plots in Figure ID (b) IRIS slit-jaw 
image in the Si IV line (~80,000 K). The image is processed to remove the dark slit with 
partial substitution of the previous-time image, (c) IRIS slit-jaw image in the Mg II line 
(10,000-15,000 K). The same process was applied as for the Si IV image, (d) SDO/AlA 
image in the Fe IX passband image (~1 MK). An active region (AR 11877) is located on 
the right-hand side out of the held of view (FOV). Bright loops that are superposed in front 
of the prominence along the LOS come from the edge of the active region. See the online 
movie A for time series of these images. (An animation of this hgure is available.) 
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Fig. 2.— Time variation of LOS velocity and line width and height dependence of oscil¬ 
lations. (a) ifmode/SOT image in the Ca II H line. The solid white line is the location of 
the IRIS slit (the second slit from the top shown in Figured]), (b) Space-time plot of LOS 
velocity derived from Mg II k spectra along the slit shown on panel (a). The notations L and 
H on the left axis correspond to the height locations indicated by the same ones on panel (a), 
(c) The same as panel (b), but of line width. Note that the values on panel (c) are the full 
width at half-maximum (FWHM). (d-e) Height dependence of LOS velocity and line width 
binned for different height ranges separated by 10" (7,200 km on the Sun) indicated on panel 
(a). The dashed and dotted lines on panel (e) indicate the instrumental line broadening (5.5 
km s“^ in FWHM) and combination of the instrumental and thermal broadenings (6.9 km 
s“^ in FWHM), respectively. 
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Fig. 3.— Heating of prominence threads. Left: Space-time plot of Hinode/SOT Ca II 
image (~10,000 K) along the green box shown in Figure [H The intensity in the green box is 
summed vertically to include the entire part of oscillating threads. Right: The same as the 
left panel, but for IRIS Si IV image (~80,000 K). Several paths (indicated by circles) show 
the transition of temperature from cool to hot, as the threads in the Ca II line fade away 
and co-spatial threads appear in the Si IV line. 
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Fig. 4.— Examining phase difference between transverse displacement and LOS velocities. 
Space-time plots of iLmode/SOT Ca images at four locations shown by the red bars (P1-P4) 
in Figure [U The yellow crosses are the central positions of the threads. The red diamonds 
indicate the LOS velocity derived from Mg II k spectra at the corresponding positions of 
the thread center, (a) PI shows clear oscillations with 180° phase difference between the 
transverse motions and LOS velocities, (b) P2 presents similar, but with 90° difference, 
(c-d) P3 and P4 show phase difference between 90° and 180°. 
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Fig. 5.— Phase difference between transverse displacement and LOS velocities at neigh¬ 
boring slit locations. Similar to Figure HJ but for all LOS velocity information at and near 
the oscillating threads. The slit locations are shown by colored lines on the upper panels. 
The colored diamonds are the LOS velocities along the respective colored slits. We can find 
some coherence of the velocity patterns over a significant distance away from the threads. 
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Fig. 6.— Synthetic intensity and LOS velocity from a simnlated prominence flnx tnbe with 
a transverse oscillation for a LOS angle of 45° in the plane perpendicnlar to the tnbe axis, 
(a) Intensity of the Mg II k line with a spatial resolntion of O.li? {R being the radins of 
the tnbe). (b) Intensity of the Si IV line with a spatial resolntion of 0.3R. (c) Space-time 
plot of the LOS velocity for a slit perpendicnlar to the tnbe axis located at the center of the 
tnbe with a spatial resolntion of 0.3R. An arrow shaped pattern develops dne to resonant 
absorption, (d) Doppler signals along the green, orange and red slit locations, with respective 
colors (whose colored edges are shown in panels (a) and (c)), centered at a distance of 1.5i?, 
OR and —1.5R from the tnbe’s center, respectively. The black dashed line corresponds to 
the transverse displacement of the tnbe in the POS, calcnlated from Ganssian fits to the 
intensity image (a) for each time step. Times withont Doppler signal correspond to regions 
for which the signal is too low. The phase differences between the LOS velocities and POS 
displacements match those in the observations (FignreH]). 
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Fig. 7.— Schematic representation of the physical model explaining our observations and 
other considered interpretations, (a-b) The transverse MHD wave model: a single flux tube 
and multiple flux tubes subject to a kink wave, (c-d) Flux tube with a torsional Alfv?n 
wave alone or coupled with a kink wave, (e-f) Flux tube with rotation alone or coupled 
with a kink wave. The observed threads can each correspond to a flux tube (a,d,f) or to 
high density regions within a flux tube (b,c,e). Locations of the IRIS slit are sketched as 
green, yellow and red boxes, as in the observations ([5]). The width of the slit corresponds to 
that of the thread. The predicted Doppler signal at each slit location (solid) and the POS 
motion (dashed) in time is sketched on the right side in each panel. For (c), two cases can 
be distinguished depending on the location of the slits over the flux tube. In each case the 
middle position is the equilibrium position for an oscillation. In the top/bottom position, 
the thread is the furthest away/closest from/to the observer. For the transverse MHD wave 
model (a,b) the azimuthal motions result in a blue/red shift at the topmost/bottom position 
(blue/red arrows), as observed. None of the other scenarios match the observations. 
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Fig. 8.— Source function and optical thickness in Mg II k and Ca II H lines. Source 
function in Mg II k 2796.35 A (left) and Ca II H line center (right) for a cross-section of 
the tube along the tube’s center. The yellow curves correspond to contour lines for different 
values of the optical thickness (stated in the legends) as seen from the top (as indicated by 
the arrow), corresponding to an LOS angle of 0°. The axes are set to match those in the 
numerical simulation in Paper 2. 
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Fig. 9.— Summary of all the physical scenarios considered. Stating the expected LOS 
velocity phase with time (together with the changes of sign, denoted as plus or minus sign) 
at the two green and red slit locations considered in the previous sketches, and the expected 
phase difference between the POS motion of the thread and the LOS velocity, x/y nota¬ 
tion indicates phase differences of x to y. The hrst column denotes the observations. Our 
numerical modeling applies to the last two columns. These are the only cases among those 
considered above which correctly match the observed features. 














